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Abstract 

We propose a method to characterize discrete time evolution equations, which gener- 
ahze discrete time soHton equations, including the g-difference Painleve IV equations 
discussed recently by Kajiwara, Noumi and Yamada[|l[|. 

1 Introduction 

A discretization of independent variables of integrable nonlinear differential equations 
breakes their integrability in general. Therefore it is remarkable that there exist certain 
discrete analogue of integrable differential equations which preserve integrability [p| ~ ||T6| . 



Integrability of ODEs can be tested by studying whether exist singularities which depend 
on initial values, a method called Painleve test. Therefore it is natural to look for a discrete 
analogue of the Painleve property which enables one to predict behaviour, deterministic or 
chaotic, of a sequence of map starting from some initial value. There has been, however, 
not known a discrete analogue of Painleve test despite of some useful proposals [0]. 
In our previous paper we have studied a discrete analogue of the Lotka-Volterra equa- 



tion which is known completely integrable |T^. Considering it as a sequence of a map it 



was observed that at every step of the map there are two possibilities to be chosen. Hence 
the orbit is not deterministic. Nevertheless the map is compatible with integrability of the 
equation. This happens due to the fact that one of two types of the map, which we call 
B-tj^e, does not generate a new orbit but simply exchanges two different orbits created 
by the other type (A-type) of the map. 

Starting from a symmetric discrete Lotka-Volterra equation (dLV) under the periodic 
boundary conditions the map is determined solving a quadratic equation. The expression 
under the square root which appears in the solutions turns out to be a perfect square in 
this problem, hence no singularity arises in the map. Out of the two maps the A-type 
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map becomes one which reduces to the ordinary Lotka-Volterra equation in the continuous 
time hmit, whereas the B-type map has no correspondence in the continuous time. If the 
equation is deformed arbitrary, the square root singularity remaines and one can not 
obtain regular map. 

We now ask if the disappearance of the square root singularity is special in the dLV. 
The purpose of this paper is to show that the phenomena is quite general among various 
discrete integrable systems including discrete KdV, discrete time Toda and discrete KP 
equations as well as their g-difference versions. 

In order to clarify our problem we review briefly, in the following section, the argument 
of the dLV in our previous paper We generalize the dLV in §3 and show that the 
same argument holds in this case too. This generalization includes, as a special case, 
the g-difference fourth Painleve equation discussed recently in ref[|l|. We can further 
generalize the scheme including arbitrary functions which we discuss in §4. g-difference 
discrete KdV, Toda, KP equations will be obtained as special reductions. In the last 
section conservation laws are discussed. 



2 Discrete Lotka-Volterra equation 

The discrete time Lotka-Volterra system is defined by the set of equations 

xi^\l-5xt\)=xi{l-5xi^,), jeSZ, n = l,2,---,N (1) 

together with the boundary conditions 

x^+i = x\, j e6Z. 

Here 6 denotes the minimum interval of the time step. This system is known completely 
integrable[Q. The integrability is justified by noting that it can be obtained from the 
discrete time Toda lattice via a Miura transformation, or existence of sufficient number 
of conserved quantities. 

We have investigated behaviour of solutions to this equation and found two types of 
map, one of which has a correspondence in the continuous time limit and the other does 
not|jT8|. To be more precise let us consider the case of = 3, and write the equations as 

(2) 
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-^) 


= X 
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-y) 


' y 


;i 


— x) 


= y 


[1 


-^) 




\ 




= z{ 







where we use the notations 

(a;, z) = {5x{, 6xi, 6xi), (x, y, z) = {5x{^^, 6xi'^\ Sxi'^^). (3) 
In order to clarify our problem let us consider the following equations at first 

x{l-z)=X, y{l-x) = Y, z{l-y) = Z. (4) 
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Here (X, Y, Z) are some functions of (x, z). Solving for x, we find 

I + X -Y - Z ± J {I - X - Y - Zf - AXYZ 



X 



(5) 



2(1 -Z) 

and similarly for y and z. 

This result shows that the mapping is not single valued. Hence there arises naturally 
a question, how this double valuedness is compatible with the integrable nature of the 
Lotka-Volterra equation (^. This problem is answered if we substitute the right hand 
sides of (H) into (X, Y, Z) of @. In fact we obtain 



I X \ 







1 — z + xz 






1 — z + zx 


fx\ 




1 — X + yx 




1 — X + xy 


y \ = ' 




z 

\ 1-y + zy 






,1-xj 



A - type 



(6) 



: B — type. 



Note that the B-type map does not generate orbits but simply exchanges the variables. 
If we operate B six times, for instance, all points return to the original place after the 
sequence of the map: 



{x,y,z) 



z,l 



x) -> {z,x,y) (1 

(y, z,x) -> {1 - z,l - x,l - y) 
There exist two constants of the equation of motion (^: 

C 
D 



x,l-y,l 
{x,y,z). 



X + y + z — xy — yz — zx 
xyz{l - x){l - y){l - z). 



(7) 



They are conserved by both A- and B-type of maps. Moreover we also found that there 
are two constants 

r = xyz, s = {1 — x){l — y){l — z) (9) 
which are conserved under the A-type map but are exchanged by the B-type map: 



(r, s) 4^ (r, s) 
(r, s) 4^ {s,r) 



for A — type, 
for B — type. 



In order to see the continuous time limit of the map we must recover S in the old 
variables (|^). We will find that in the limits 6—^0, j t, the equations (0) as well as 
the A-type map of (||) reduce to the symmetric Lotka-Volterra equations 



dx ^ ^ dy 
dt ' dt 



dz . 
-dt='^y 



y{^-z)^ — = z{y-x)- (10) 

On the other hand the B-type map does not have a proper limit as it was discussed in 

m. 
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3 (/-Lotka-Volterra maps 

The integrability of the Lotka-Volterra map (Q) owes to its special symmetry under the 
cychc permutation of variables. If we replace it by, for instance, 

x{l — z) = ax{l — y) 

y{l - x) = az{l - z) (11) 
. ^(1 -y) = az{l - x), 
the integrability is lost when a ^ 1. Therefore it is highly nontrivial question whether 
there exists generalization of the discrete Lotka-Volterra system @j which preserves inte- 
grability. 

We consider in this section the following set of equations as a generalization of (Q): 

x{l — z/c) = aX, 

y{l-x/a) = hY, (12) 

. z{i - y/h) = cZ, 

where (X, Y, Z) are some functions of (x, y, z) and a, 6, c are arbitrary constants. By the 
same reason we discussed in the previous section, we expect maps of the form (^) when 
we solve (^) for (x, |/, z) in general. In fact we obtain 



X 



l + X-Y -Z±J{l-X-Y -Zf- AXYZ 



(13) 



2(1 -Z) 

Motivated by our previous experience we try to find (X, Y, Z) such that the expression 
under the square root is a perfect square. We found the following three cases satisfy this 
condition: 



■ x/a) 

y/b) 

z/c): 

x/a) 

y/b)- 

-Z/C): 

x/a) 

y/b)- 

Solving them for (x, |/, z) we obtain 

I -by 



I 2 

a X- 



a^x(l — by), 
-- Py{l - cz), 
c^z{l — ax), 

abyil — cz), 
■ hcz{l — ax), 
cax(l — by), 

acz{l — ax), 
: 6ax(l — by), 
cby{l — cz). 

bcyz \ 



(14) 
(15) 
(16) 



b'y 



1 ~ cz + acxz 
1 — cz + acxz 



1 — ax + abxy 
2 1 — ax + abxy 



c z 



: A — type 



V 1 — 6y + bcyz I 



(17) 



/a(l 


-by) 


6(1 


— cz) 


\c{\ 


— ax) 



: B - type 



for the first set, 







^ o-by 






bcz- 








y \ = ' 











1 — cz + caxz ^ 

1 — ax + abxy 
1 — ax + abxy 

1 — by + bcyz 
1 — by + bcyz 

1 — cz + acxz 



A - type 



(18) 



a(l — cz) 
6(1 — ax) 
c(l - by) 



: B — type 



for the second set, and 



/ 1 — ax + abxy \ 

aCZ- ; ; 

1 — by -\- bcyz 
1 — by + bcyz 



abx 



bey 



I — cz + acxz 
1 — cz + acxz 



A — type 



\ 1 — ax + abxy I 

: B — type 



(19) 



a(l — ax) 
b{\-by) 
c(l — cz) 



for the last set. 

Note that the first set (0) is the Lotka-Voherra map (^ when a = b = c = 1. The 
A-type map (0) of the second set is the g-Painleve IV studied by Kajiwara, Noumi and 
Yamada[|l[. 

From these expressions it is apparent that the above three cases are symmetric under 
the cychc permutations: 



(ax, by, cz) — ^ {by, cz, ax) — > {cz, ax, by). 



(20) 



We will call this symmetry the cyclic symmetry. Owing to this symmetry it is sufficient 
to consider only one of three equations. 

In addition to the cyclic symmetry, the symmetry under the exchange of the variables 

{x/a,y/b,z/c) {ax,by,cz) 

in equations (|I1|),(|TB|),([TB|) enables us to solve them for {x,y,z) in terms of {x,y,z). We 
call this symmetry the time reversal symmetry. For instance we solve, corresponding to 

(0), 



a'^x{l — cz) = x{l — y/b), 
b^y{l - ax) = y{l - z/c), 
^ c^z{l — by) = z{l — x/a). 



(21) 
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for (x, ?/, z) and obtain 



f X l — z/c + yz/bc\ 

a? 1 — y/h + xy/ab 
y 1 — x/a + xz/ac 

z 1 



■ z/c + yz/bc 
y/b + xy/ab 

\ 1 — X / a + xz / ac J 



A - type 



ni-z/c)/a\ 

{l-x/a)/b : B - type. 

\{l-y/b)/c] 

We could obtain this result by solving (p!7|) for (a;, |/, ^). 

Generalization of our argument to > 3 is straightforward. For example the 4-point 
discrete Lotka-Volterra equation is generalized to 

' x{l — w/d) = a'^x{l — by) 
y{l — x/a) = b'^yil — cz) 
z{l - y/b) = c^z{l - dw) ' 
. w{l — z/c) = d'^w{l — ax) 

and is solved by 



fx\ 

y 

z 



1 - 


- by - 


- cz + bcyz + cdzw 


1 — 

1 - 


cz - 
cz — 


dw + adxw + cdzw 
dw + adxw + cdzw 


^- 

1 - 


ax - 
ax - 


- dw + abxy + adxw 

- dw + abxy + adxw 


'l- 

2 1- 


- ax 

- ax 


-by + abxy + bcyz 
— by + abxy + bcyz 


1 - 


- by 


- cz + bcyz + cdzw 



V 

/ a(l - by) \ 
b{l - cz) 
c(l — dw) 
. V d{l — ax) J 



: B - type. 



(22) 



A — type 



(23) 



4 Further generalization 

Let us further generalize (O) to 



' x(W - z/c) = aX, 
< y(U — x/a) = bY, 
j{V-y/b) = cZ, 

where U, V, W are some functions of x, y, z. (|13D is then changed to 



(24) 



UVW - WY -UZ + VX± JiUVW - WY -UZ - VXy - AXYZ 
^ = " 2iVW - Z) • ^2^) 
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We want to know X, Y, Z such that the expression under the square root in ( P^D becomes 
a perfect square. We find that 

' x{W - z/c) = a?x{U - by), 

< y{U -x/a) = b'^y{V -cz), (26) 
jlv-y/b) =c^z{W-ax), 



and two other equations, obtained from (^) by cychc permutations (pO]), do satisfy the 
requirement. Solving (|26| ) for {x,y, z) we obtain two maps 



fxX 

y 



/ , UV - bVy + bcyz \ 

(XX 

VW - cWz + acxz 
, VW - cWz + acxz 



WU — aUx + abxy 
2 WU — all X + abxy 



\^ ^ UV -bVy + bcyz J 



: A — type 



(27) 



a{U - by) 
b{V - cz) 
c{W - ax) 



: B - type. 



The maps of ( p7[ ) have constants 



R : = 



xyz 



abcxyz ' 
under the A-type map and 



S :-- 



{aU -x){bV -y){cW - z) 
abcxyz 



{U -by){V -cz){W -ax) 
{aU - x){bV - y){cW - z) 



xyz 

{U -by){V -cz){W -ax) 



[28] 



under the B-type map. Note that the B-type map exchanges R ^ S. 

When {U,V,W) are functions of {x,y,z), the time reversal symmetry is lost, hence 
[x, y, z) obtained by solving (^) will not be rational polynomials of {x, y, z) in general. 
The symmetry is recovered if U, V, W have some symmetry under the exchange {x, y, z) ^ 
[x, y, z). Apparently the g-difference Lotka-Volterra system satisfies this condition. We'll 
present other examples in what follows. 



4.1 g-KdV maps 

We have noted that the A-type map of 
U, V, W are given by 

U = V = W = abcxyz 
the requirement is fulfilled. Namely we can substitute another expression 



preserves R, irrespective of U,V,W. Thus if 

(29) 



U = V = W = xyz/abc 
of U, V, W to the left hand side of (|26|) to solve the inverse map. 



(30) 
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In fact a substitution of 



into the A-type map of ( JI7\ ) yields 
abxy + a?hcx^yz \ 



bcz 



cax- 



1 — bcyz + ab'^cxy'^z 
1 — bcyz + ab'^cxy'^z 

1 — cazx + abc^xyz"^ 
1 — cazx + abc^xyz'^ 



(31) 



1 — abxy + a%cx'^yz ) 
We can solve this map inversely for (x, 2) and obtain 



z 


a%c 


— abxz + \ 


a2 


ab(? 


— acyz + x|/^2 


X 


ab^c 


— bcxy + xy^^ 




d^bc 


— abxz + x^l/z 


V 


ab(? 


— acyz + xyz^ 


C2 


ah^c 


— 6ca;?/ + xy'^z / 



It is not difficult to convince ourselves that the same map is obtained by solving 
(x, y, z) after the substitution of (|30D. 



for 



Now we like to show a relation of (0) with the discrete KdV map. The latter has 
been proposed in the form]^, |T^ 



(32) 





1 




1 


X — 


y 


= X — 


z 




1 




1 


< y- 


z 


= y- 


X 




1 




1 


z — 


X 




y 



A direct substitution of 



/ 1 

y 



xy + x'^yz \ 



1 — yz + xy^z 
1 — yz + xy^z 

1 
1 



X- 



zx + xyz^ 
zx + xyz^ 



(33) 



will show that this solves 
From this reason we call 



21). But (13 



V 1 — xy + x'^yz / 
is simply a special case of (|3T|) with a = b 



X b 1 

= ax , 

ay cz 

- = hy , 

b z ax 

z a 1 
- - - = cz - — , 
c X by 



(34) 



the g-KdV map, which is obtained from ( ^6|) after substitutions of (|29|) and (pOl). 
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4.2 g-Toda maps 

In addition to x,y,z we introduce variables u,v,w and assume for U,V,W 
following relations: 



u 


= by + a'u 


V 


= cz + b'v 


w 


= ax + c'w 



X 

a 

y 



V 

w 
z u 

- + -■ 



+ 



Then we find solutions to (p6D 

/ 



(x\ 

y 

z 
u 

V 



, a'b'uv + ca'uz + bcyz \ 

o-^TT, r, 

b'cvw + abvx + acxz 



b'y— 



b'c'vw + ab'vx + acxz 



a'c'uw + bc'wy + abxy 
r, a'c'uw + bc'wy + abxy 

a'b'uv + a'cuz + bcyz 
b'c'vw + b'avx + acxz 



c z 



a''^U ; ; 

a'b'uv + a'cuz + bcyz 

, ,„ a'c'uw + bc'wy + abxy 

b V- ; 

b'c'vw + b'avx + acxz 

,„ a'b'uv + a'cu2; + 6cyz 

\ ^ / 

\ a'c'uw + bc'wy + a6xy / 



/ aa'u\ 
bb'v 
cc'w 
aa'x 
bb'y 

V cc'z J 



B - type. 



Note that this system is equivalent to the equations 



X 



y 



w 



: A - type 



u 



by + a'u = — h 77, cz + b'v = - -\ — -, ax + c'w = - H — - 



xu 
aa' 



aa'xu. 



yv_ 
bb' 



bb'yv, 



zw 
cc' 



cc zw, 



which is nothing but the discrete time Toda equations [0 when a = b = c 
c' = 1. Hence we obtained a q difference discrete time Toda system. 



4.3 g-KP maps 

We now consider three sets of equations of the type (^Bf) : 



ymiUm •^m/ O'nri) ^m2/"i(^^ Cm,Z-> 
, Zm(^m ym / bm) Cjj^Zfji(\Vn 



Om-^m ) 



m = 1,2,3 
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By introducing new variables {um,Vm,Wm) and {um-,Vm-,Wm, 



we assume for (f/, V, W) 



the following form 





+ 












Um 


+ 








bm 


d 




+ 


Um 






Cm 







TT •^m . ^m 7 . / 

^m. ~r ~r. "mVm ' '^-^^j^^m+l; 

Cm^m ~l~ bjy^_i_-^Vm+l} 

- (Xm^m ~l~ Cjy^j^-^Wm+l- 

Here the periodic boundary conditions Xm+j, = Xm,Um+3 



m 



1,2,3 



(39) 



Ur 



ttm+s = CLm etc. are used. 



Solving (BHI) for 



Wm) we find 



^ Xm 


\ 


Um 
















\Wm 


) 



I 



2 



Q'm+l^m+l^'rn+l'^m+l ~l~ (^rn+l^rn'^m+l^m ~l~ bmC^mUm^m ^ 
^m+l'-'m+l^"i+l^m+l ~l~ ^m+l'-^™^m+l'^'Ti (^mCm^m^m 



2 ^ ^m+l'-'m+l^™+l^'Ti+l ~^ ^m+l'^™^"^+l'^'" Cm^m^mXr, 



~r C^j^i^Om'U^ m+lUm ' "m^m-^mUn 
C^j^l^Cl^j^lWm+lUm+l ~l~ (^m+l^m'^m+iym ~l~ (^mbm-^mUn 



a' 



m+l^m+l^m.+l^m+l (^m+l^ni^m+l^m ~l~ ^m^^mym^m 

I I ^m,+l'-'m+l^''i+l^"i+l ~^ ^m+l'^™^"i+l'^"^ ~^ 0,m(^m-^m^', 

^m+l^m'^m+1 ~ 7; ~ ; ~T 

^m+l^m+l^m+V^m+l ~r dm+l^m'^'m+l^m ~r OmCmym^', 

,/ C^_j_-|^(3j^^-|^'UJ,^-j_l1tn^-f 1 ~l~ C^_^_^bm'Wm+iym ~l~ (^mbm-^my 



m 
m 



^m+l^m'^m+1 , 



^m+l'-'m+l^'Ti+l'U^m+l ~l~ ^m+l'-^m^m+l'^'Ti ~l~ Cm(^m^mXm 
(^m+lbm+l^'m-+^^fTT-+^ ~^ '-^m+lCm^m+l^m ~l~ bmCmym^m 



(40) 



-I'-^m+l'^m+l^m+l ~l" Cjji_^_ibmWm+iym (^mbm^mym / 

■ Am -type 



^m^m'^n^ 



bmbfYi^ym 



V 



Bm-type 



with m = 1,2, 3. All eight combinations of the A^-type and B^-type maps 

(Ai, ^2, As), (Ai, ^2, 53), • ■ ■ , (5i, 52, 53) 

are possible map. 

We would like to note here that the set of equations ( pHf ) together with ( pOf ) are 
g-difference discrete time KP-equation or the discrete 2D Toda equation. In fact the 
substitution of (P^) into ( PB]) shows that the set of equations (^Sp and (^^ are equivalent 

to 

(41) 
(42) 



X7- 



+ 



XmnUmn 
(^mnbmn 
'^m,n+l 



(^m,nbm+l,nXm,n'^m+l,n^ 
m,n+l-^m,n+l ~l~ ^m+l,n^»n+l,n) 
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respectively. Here we use the notations 

(•^771,1; '^m, 2? -^m, 3) {p^miymi^m)i (^m,l) ^m,25 ^771,3) i^'^m^ '^mi ^rri) 1 
(^0'm,li 0"in,2i 0"m,3^ ('3'm; ; Cm) > (^m,l) ^m,2) ^01,3) ('^m) ^m' '^m) ■ 

The discrete KP (or discrete 2D Toda) equation has been given by|lTO| 

+ + + 1) = ^) + v)- 

The correspondence will be established if we put 

and (^) and identify the variables according to 

u = m + n, = j + fn, rj = —m. 

Therefore (^) together with (|39D is a g-differnce discrete KP equation. 

5 Discussions 

We have found a large class of maps which share the same properties with those of discrete 
Lotka-Vorterra map. All of them have two types of map, which we called A-type and B- 
type. The B-type map does not create a new orbit, but exchanges orbits generated by 
the A-type map. 

In the rest of this paper we'll discuss about conservation laws in our system. We 
already pointed out in the case of dLV that C and D of and (|) are conserved along 
the sequence of the map. If we restrict a map to the A-type there are additional constants 
r and s given by @. 

This is true also in other cases ii a = h = c = 1 and U = V = W . In fact from (|26|) 
we see that 

Cu = {x + y + z)U — xy — yz — zx, Djj = xyz{U — x){U — y){U — z) 

are constants. The dLV and dKdV belong to this case. If we consider only A-type map 
and assume a = h = c = 1, (pT]) shows that r = xyz is a constant of the map irrespective 
to U, V, W. 

The discrete time Toda lattice contains more variables. We can read off conserved 
quantities from (|37|) as 

u + v + w + x + y + z, ux, vy, wz 

ii a = h = c = a' = h' = d = 1. In the case of A-type of map 

xyz, uvw 
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in (51 



are also conserved as long as a = b = c = a' = b' = c' = 1. 

When the maps are generalized by introducing the parameters a, 6, c, • ■ •, these quanti- 
ties are not conserved along the map anymore. For instance we find, in the case of A-type 
map of (p?!), r = xyz changes according to 

xyz = q'^xyz, q = abc. (43) 

This means that the scale is changed by an equal rate at every step of the map. Hence 
the quantity xyz is not constant. 

The scale change of this form was identified with a time evolution of the system in 
ref[|[]. From this point of view we can interprete 2 In g as a velocity of the variable \n{xyz). 
Namely, writing 

r{t) = xyz = g^V(O), 

we have 

d\nr(t) 

As we already know this is a constant of the map of A-type. In the special case of g = 1, 
the velocity is zero, hence r{t) itself turns out to be a constant. Therefore it is natural 
to understand that the quantity R in (|28|) is a conserved quantity of the A-type equation 
of motion of (pT]). ln(gi?) is a difference of variables at two subsequent steps of the map. 
One may interprete it as a center of mass velocity of the coordinates (Inx, Iny, Inz): 

ln(gi?) = (In a; -|- In 1/ -|- In ^) — (In a; -|- In ?/ -|- In 2;) = —{lnx + \io.y + In z) 

In addition to R we have another quantity S which remains constant under the map. 
S, which is given by (|28|) , also depends on variables of two steps of the map. It has no 
direct interpretation like R, but becomes R by the B-type map. 
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